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Q ■ Abstract 

The theory of resonant Auger decay of atoms in a high intensity coherent X-ray pulse is presented. 
The theory includes the coupling between the ground state and the resonance due to an intense 
X-ray pulse, taking into account the decay of the resonance and the direct photoionization of the 
ground state, both populating the final ionic states coherently. The theory also considers the 
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, impact of the direct photoionization of the resonance state itself which typically populates highly- 

O ' excited ionic states. The combined action of the resonant decay and of the direct ionization of the 

C/3 ■ 

ground state in the field induces a non-hermitian time-dependent coupling between the ground and 
the 'dressed' resonance stats. The impact of these competing processes on the total electron yield 
and on the 2s^2p^(^D)3p spectator and 2s^2p^ participator Auger decay spectra of the Ne 
' Is— T-Sp resonance is investigated. The role of the direct photoionization of the ground state and 
■ of the resonance increases dramatically with the field intensity. This results in strong interference 

! effects with distinct patterns in the electron spectra, different for the participator and spectator 

I final states. 
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I. INTRODUCTION 



The resonant Auger (RA) effect of atoms, discovered more than 30 years ago [l|, |2|], has 
been extensively studied over the past decades experimentally by utilizing conventional syn- 
chrotron radiation sources and theoretically (see, e.g., the review [sl] and references therein). 
In contrast to X-ray photoelectron spectroscopy (XPS), the RA decay spectra are recorded at 
the exciting-photon energy tuned to a resonance. In the vicinity of core-excitations of atoms, 
the resonant photoionization channel (i.e., the excitation and Auger decay of a highly-excited 
electronic state) dominates over the direct channel (i.e., non- resonant photoionization) pop- 
ulating the same final ionic state, owing to the large energy of the exciting X-ray photons. 
As a result, the RA decay electron spectra provide uniquely important information on the 
electronic structure and decay mechanisms of highly-excited electronic states. 

The impact of the direct photoionization on the RA spectra is nevertheless not negligible. 
There are several theoretical and experimental studies indicating the relevance of the weak 
direct photoionization channel for the interpretation of RA spectra of atoms and molecules 
4-16i. Perceptible fingerprints of the interference between the dominant resonant and the 
weak direct photoionization channels were observed in the angular-averaged RA decay spec- 
;ra and in the branching ratios of different decay channels scanned across the resonance 
^-|9|. This interference manifests itself even more distinctly in the angular-resolved RA 
decay spectra I^IG]. As demonstrated in these works, the weak direct photoionization 



channel leads to a broad photon-energy dependent dispersion of the RA electron angular 
distribution and residual ion polarization parameters. 

The advent of X-ray Free Electron Lasers (XFELs) HQ raises the fundamental ques- 
tion of how the well-studied processes of interaction of matter with electromagnetic field will 
be modified under extreme field conditions. Compared to conventional synchrotron radia- 
tion sources, the intensities of the electromagnetic fields generated by XFELs are by several 
orders of magnitude stronger, and the durations of the light pulses can be made to be similar 
to the typical lifetimes of highly-excited electronic states (typically a few femtoseconds). The 
first theoretical studies of RA decay of atoms exposed to strong X-ray pulses |l9l-l2l| showed 
that the stimulated emission from the resonance back to the ground state starts to compete 
with the Auger decay. The interplay between the resonant excitation and stimulated emis- 
sion results in Rabi oscillations between the ground state and the resonance within its Auger 
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decay lifetime and pulse duration [l9|, l20|, which leads to spectacular modifications of the 



RA spectra. It has been demonstrated in Refs. [20|, |2l| that in the presence of such strong 
electromagnetic fields the direct photoionization of the ground state may start to play an 
important role. 

Obviously, a final ionic state is populated coherently by both the direct photoionization 



and the resonant Auger decay thus natural 
spectra. On the other hand, as discussed in 



y inducing interference effects in the electron 
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2l| , there will be leakage of the population 



of the ground state by the direct photoionization into all possible final ionic states (total 
photoionization). We shall demonstrate here that during the duration of a strong pulse the 
direct photoionization of the resonant state itself to produce preferably highly-excited ionic 
states may also take place and compete with its decay via Auger and stimulated emission to 
the ground state. It has not yet been shown and understood how all these processes evoked 
by a strong field infiuence the total electron yield and the spectra and filling this gap is the 
major goal of this paper. 

The derivation of the theory of electron spectra of atoms exposed to strong X-ray pulses 
is presented in Sees. Ill AHII Dl The readers not interested in the details of the derivation are 
referred to Sec. Ill Et where the final set of equations is collected and analyzed. In Sec. IIIIl 
the derived theory is applied to the Auger decay of the ls^2s^2p^3p^ resonance of the Ne 
atom into the conceptually different spectator and participator decay channels. We conclude 
with a brief summary. 

II. THEORY 

The currently operating XFEL, Linac Coherent Light Source (LCLS), does not produce 
so far a monochromatic radiation and its X-ray pulses consist of many spikes with random 
fluctuations of the frequency, phase and amphtude |171] . The impact of these problems on the 
RA effect has been studied in In addition, the X-ray pulse properties are compressed 



through the propagation in a resonant medium |2l[ |. 

In the present work we concentrate on the physics a single atom undergoes when exposed 
to a coherent and monochromatic X-ray pulse. The knowledge of this physics is prerequisite 
for further studies. For simplicity we assume a pulse linearly polarized along the z axis with 
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FIG. 1: Schematic representation of the resonant Auger effect of an atom in weak (left panel) and 
strong (right panel) fields. Left panel: In a weak electromagnetic field, the resonance is excited 
and decays populating a final ionic state which is also populated by the direct photoionization of 
the ground state. Both pathways coherently superimpose, but due to the high X-ray frequency, 
the resonant pathway dominates by far. After the pulse expires, the residual population of the 
ground state remains nearly 1. Right panel: In the case of a strong field, the ground state and the 
resonance are strongly coupled by the field. The usual direct coupling which is known to induce 
Rabi oscillations between these two states is modified by an interesting term which is indicated by 
a vertical helix double-arrow. The ground state, the resonance and also the produced ionic states 
are subject to time-dependent losses (leakages) due to photoionization as indicated by the inclined 
helix arrow from the respective state. The modified coupling and the leakages influence the RA 
effect and lead to a damping of the Rabi oscillations. After the pulse expires, the population of 
the ground state can be very small. 



the field: 

£{t) — £o{t) cosLot — So g{t) cosut. (1) 

Here So is the peak amplitude, and the pulse-shape function g(t) varies slowly on the 
timescale of 2t:/u}. The cycle-averaged intensity of the field is given in atomic units via 
(1 a.u. = 6.43641 xlO^^ W/cm^) 



(2) 



where a — 1/137.036 is the fine structure constant. 



The processes relevant for the present study are schematically shown in Fig. [T] In a 
weak electromagnetic field (left picture), the resonance is excited and decays populating a 
final ionic state which is also populated by the direct photoionization of the ground state. 
Both pathways coherently superimpose, but due to the high X-ray frequency, the resonant 
pathway dominates by far. As the pulse expires, the non-resonant pathway expires as well, 
and the population of the ground state remains close to the unity. In the case of a strong 
field (right picture), the ground state and the resonance are strongly coupled by the field. 
We shall demonstrate below that the usual direct coupling which is known to induce Rabi 
oscillations between these two states is modified by an interesting term derived in Sec. Ill CI 
and discussed in Sec. Ill El The modification of the coupling between the states is indicated 
in Fig. [T] by a vertical helix double-arrow. In addition, there are two mechanisms of losses 
relevant to the RA effect in strong fields: The total photoionization of the ground state and 
of the resonance lead to a leakages of the population of the respective states, as indicated in 
Fig. [T]by the inclined helix arrows from the ground state and from the resonance, respectively. 
These two loss mechanisms have interesting impact on the RA effect. The leakages and the 
modified coupling all lead to damped Rabi oscillations of the population between the ground 
state and the resonance. Of course, these processes, except of the Auger decay itself, are 
operative only when the pulse is on. After the pulse expires, the final population of the 
ground state can be much smaller than 1. 

In the presence of a strong pulse, multiple ionization of the atom takes place as well 



18| . There are several mechanisms responsible for that. One is the ionization of the final 



ionic states populated via the RA effect (indicated by the inclined helix arrow from the ionic 
state in Fig. [1]). The ionization of the resonance mentioned above is another. This ionization 
produces preferentially highly-excited ionic states, most of them with a core electron missing 
as in the resonance state. These ionic states are then likely to undergo Auger decay producing 
thereby doubly-ionized or even higher ionized atoms. All the processes discussed here can, 
in principle, be experimentally separated from each other by measuring the energies of the 
ejected electrons. 

Below we concentrate on the processes shown in the scheme of Fig. [H We mention that 
the leakage from the ionic states is not relevant for the RA effect as long as one measures 
electrons and does not intend to measure the ions. To describe the RA effect theoretically 
as a function of time, one needs to solve the time-dependent Schrodinger equation for the 
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atom and its interaction with the field (atomic unis e = rUe = h = 1 are used throughout) 

z^(t) = Hm{t) = (#0 + zS{t)^ ^{t). (3) 

rn 

The present approach to solve Eq. is similar to that reported in [22]. For transparency 
of presentation we derive below step by step the equations for the different mechanisms 
depicted in Fig. [1] and collect all results and discuss them in Sec. IIIEI 



A. Resonant channel 



As a first approximation, we include the resonant channel, consisting of the excitation 
and Auger decay of the resonance and the stimulated emission from the resonance. We start 
with the following ansatz for the total wave function of the three-energy-level system: the 
ground state |/), the resonance \R), and the final ionic state plus Auger electron \Fe), with 
the energies Ej, E^, Ep and e, respectively 

^{t) = ai{t)\I) + aR{t)\R) + j apit, e)\Fe)de. (4) 

Here a/(t), aji{t), and aF(t,e) are the time- dependent amplitudes for the population of 
the I J), \R), and \Fe) levels, respectively. The transitions between the atomic states are 
described by the following matrix elements of the total Hamiltonian ([3]) 

{Rmm = (i?|5|/)^^e--* = Dnit)e~^-\ (5a) 

{Fe\H{t)\R) = V. (5b) 

In Eq. ( l5al) . the rotating wave approximation j23] has already been utilized, and, in contrast 
to the rapidly oscillating factor e"""^*, the function Dii{t) varies slowly on the timescale of 
l-n jbj. Here and below, all transition matrix elements are assumed to vary slowly with the 
energy across the resonance and are replaced by their mean values. 

By substituting ansatz (jlj) in the time-dependent Schrodinger equation and projecting 
the result onto each state, one obtains the following set of differential equations for the 
time-dependent amplitudes 

m/(t) = Eiaiit) + L>]j(t) e+^'^*Si?(t) (6a) 
%aR{t) = Dnit) e-'^*ai{t) + ERaR{t) + j V^apit, e)de (6b) 



iapit, e) = VaR{t) + {Ep + e) apit, e). (6c) 

It is possible to eliminate the rapidly oscillating factors e''^*'^* from the system of equations IQ. 
For this purpose we redefine the time-dependent amplitudes as follows 

a^(t) = anit) 6+^"^* and apit, e) = ap{t, e) e^^'^K (7) 



In the local approximation 2% l25||, the integral in the right part of Eq. (l6bl) is equal to 

j ap{t,e)de = -nr\V\'^aji{t), (8) 
as has been explicitly demonstrated in Ref. 26|] and will also be outlined in the next subsec- 



tion for the case of the direct ionization channel. Finally, when several Auger decay channels 
of the resonance into different final ionic states \Fjej) exist, Eq. (IHl) must be modified as 
follows 

J2 I Vjap^{t,e,)dE, = -tirJ^lVil'Mt) = -^r^a^W, (9) 
j j 
where is the total rate for the Auger decay of the resonance. 

In the above designations, the set of coupled differential equations for the time-dependent 

amplitudes reads 

idi{t) = Eiaiit) + D^^it) anit) (10a) 
ianit) = Dnit) ajit) + {Er - ^F^ - uj)aR{t) (10b) 
iixp^ (t, Sj) = Vj aR(t) + [Ep^ + Sj - u) ap^ (t, Sj). (10c) 



The system of Eqs. ffTOj) is consistent with the equations reported in 19|. It includes the 
coupling of the ground electronic state of energy Ei with the 'dressed' resonance state of 
energy Er—oj by the matrix element DR{t), as well as the leakage of the resonance population 
with total rate F^ into the 'dressed' continuum final states with the energies Ep^ -\- sj — u 
via the partial Coulomb matrix elements Vj. The leakage of the resonant population caused 
by the Auger decay mechanism is referred to hereafter as the ^ RA-leakage\ 

B. Direct channel 

In this subsection, only the direct population of the final ionic state iFe) from the ground 
state |/) via the transition matrix element 

{Fe\H{t)\I) = (F£|i|/)^^|^e-*"* = DD{t)e-'''' (11) 



is investigated. We apply the following ansatz for the total wave function of the two-level 
system: 

^(t) = ai{t)\I) + j aF{t,e)\Fe)de. (12) 

After redefining the apit, e) via Eq. ([7]), the set of differential equation for the time-dependent 
amplitudes reads 

idi{t) = Eiai{t) + j D^j^{t)aF{t,e)de (13a) 
^^(t, e) = Dnit) aj{t) + {Ep + e — u) ai?(t, e). (13b) 
The formal solution of Eq. (]13b|) is given by 



a pi 



{t, e) = -i / DD{t')ai{t') e-*(^^+^-")(*~*')rft'. 



(14) 



With the help of this expression, the integral on the right-hand side of Eq. fll3ap can be 
rewritten as 

D^^{t)aF{t,e)de = -i K{t - t')ai{t')dt' , (15) 

J J — oo 

with the non-local kernel 



K{t ~t') = j D\,{t) Doit') e-'(^^+^-")(*-*')rf5. 



(16) 



The kernel (fT6|) accounts for the effect of coupling the ground electronic state with the 
continuum of final states. In the local approximation 24, |25|, the integrations over the 
energy and time in Eqs. (ITBl) and (fTB]l can be performed analytically (see Appendix |X] for 
details). The final result reads 



J D^^{t) aFit,e)de = -z7r|L)^(t)pa/(t). 



(17) 



In the presence of several final ionic states \Fjej), Eq. 0171) must be modified as follows 



where Tph{t) is the time-dependent total probability for the direct photoionization of the 
ground state. This Tph(t) is identical with the "iphif) introduced in Ref. |20| (see discussion 
around Eq. (16) in this reference). 

In the above designations, the set of equations for the time-dependent amplitudes reads 



idiit) = {Ei--rph{t))ai{t) 



(19a) 



iapj (t, Sj) = (t) ai(t) + [Ep^ + Sj - u) ap^ (t, Sj). (19b) 

The system of Eqs. f|T9|) includes the leakage of the population of the ground state with the 
total probability for the direct photoionization Tph{t), as well as the direct population of the 
'dressed' continuum final states with the energy Ep^ +ej — uj via the partial matrix elements 
DD.j{t). The leakage of the population of the ground state owing to the non- resonant direct 
photoionization channel is referred to hereafter as the ^ GS-leakage\ The GS-leakage is time- 
dependent, taking place only during the pulse duration (see Eq. (fTTl) for the definition of 
the function Di:,{t)). 

C. Interference of resonant and direct channels 

In order to account for the interference between the resonant and direct ionization chan- 
nels, we apply the ansatz (jll) for the total wave function of the three-level system with the 
coupling matrix elements given by ([5]) and f|TT]) . The set of equation for the time- dependent 
amplitudes can be obtained as described above. After the redefinitions ([7]) it reads 



iaiit) = Eiaiit) + D\{t) aR{t) + j D^^{t) ap{t, e)de (20a) 

lanit) = Dnit) aj{t) + {Er - u) aR{t) + j a^(t, €)d€ (20b) 

iap{t, e) = Doit) aiit) + VaR{t) + {Ep + e - u) ap{t, e). (20c) 
The formal solution of Eq. f l20cp is given by 

ap{t, e) = -% r {Do(i!) ajii!) + VaRii!)} e-'^^^+'-'^^^'-^'Ut' . (21) 

J — oo 



Utilizing this expression and implying the local approximation the integrations on the right- 
hand sides of Eqs. flM and fl20bl) can be performed analytically, as shown in Appendix |A] 
for the case of the direct channel. The final results read 

j D^^it) ap{t,e)de = -in {^\DDit)\''ai{t) + D^^it)VaRit)^ , (22) 

J V^ap{t,e)de = -in {Dz?(t)y"^a/(t) + \VfaR{t)} . (23) 

In the presence of several final ionic states \Ejej), and in the above designations, the set 
of Eqs. (j20p takes the following final form 

idiit) = (^Ej - ^Tphiti^ aiit) + (^D^^it) - ^W^t)^ aR{t) (24a) 
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tdRit) = (^Dnit) - ^W{tyj aj{t) + (^En - '-Ta - i^j anit) (24b) 
iap^ {t, Ej) = Dd^ (t) ai{t) + Vj aR{t) + [Ep^ + Sj - uj) ap^ {t, ej), (24c) 
with the function W{t) defined via 

W{t) = 27^ J2^D,it)Vj^- (25) 

j 

In addition to the GS-leakage in Eq. fl24ap and the RA-leakage in Eq. (124bp . the system of 
Eqs. (!24l) incorporates the populations of the final ionic continuum state via the resonant 
and direct ionization channels coherently (the sum of two corresponding amplitudes on the 
right-hand side of Eq. fl24cl) ) . In the presence of the direct ionization of the ground state, 
the coupling matrix element between the ground and 'dressed' resonance states is modified 
by the time-dependent function fl25l) (cf Eqs. ffTOl) and flMl) ). In additional to the coupling 
caused by the field D^if) there is now the non-hermitian term —^W{t), induced by the GS- 
leakage mechanism. This term is referred below as the Leakage-Induced Complex coupling 
'LIC-coupling\ We note that the whole coupling matrix element is now non-hermitian. 

D. Direct ionization of resonance 

It is very unlikely that the energy tuned for the resonant core excitation will match 
the one required for further resonant excitation of the resonance. However, the energy of 
radiation is above the ionization threshold of the outer shells of the resonantly excited state. 
This may lead to a non-resonant ionization of the resonance state by preferably populating 
highly-excited ionic states (see Fig. [T]). In order to incorporate this mechanism into the 
theory, we extend the ansatz (jl]) by the following term 

J aF^{t,€^)\F^e^)de^ (26) 

and introduce the transition matrix element between the resonant state and the excited 
state of the ion plus electron \F^e^) with energies Ep^ and 

{F,s,\Hit)\R) = (F,£,|i|i?)^^e-^'^* = D,(t)e-*-*. (27) 

Similar to Eq. ([7]), the time-dependent amplitude ai?^(t, e*) must be redefined as 

aF,{t,e,) = aF^{t,e,)e+^'^'. (28) 

10 



Incorporation of the non-resonant ionization of the resonance results in the modification 
of Eq. (I24b|) . In the case of many direct ionization channels of the resonance into 1-^*^^*^) 
states, Eq. (124bp will contain the following additional term on the right-hand side 

Yl [ Dl^it)aF,^{t,e,^)de,^. (29) 
The set of Eqs. (12^ must also be extended by the following equation for the amplitude 

iap,. (t, e*^ = D^^ {t) aR{t) + {Ef,^ + e*^. - 2wj ap,^ (t, e^^). (30) 

One can decouple Eq. (!30l) from the whole set of equations. By substituting the formal 
solution of Eq. fl30l) into the integral fl29l) and implying the local approximation, one obtains 

j Dl^{t)aF^^{t,e,^)de,^ = ^ |Z},^(t)|2a^(t) = -lr,(t) a^(t), (31) 
j j 

Here, r*(t) is the time-dependent total probability for the direct ionization of the resonance. 

In the presence of r*(t), only Eq. (124bp from the set of Eqs. (p^ must be modified 

accordingly 

idnit) = (^Dnit) - '^W{t)^ ajit) + (^Er - ^[Fa + r,(t)] - u?j aR{t). (32) 

The term —^T^:(t) aji(t) on the right-hand side of Eq. (!32|) describes leakage of the popu- 
lation of the resonance due to direct ionization, referred to hereafter as ^ RD-leakage\ The 
RD- leakage is time-dependent. It competes with the two alternative mechanisms for depop- 
ulation of the resonance, the time-independent RA-leakage (due to Auger decay) and the 
time-dependent couphng with the ground state (responsible for the Rabi oscillations). 



E. Final equations and discussion 



Let us introduce the time-dependent vector of the amplitudes of the total wavefunction 
describing the population of the various involved atomic levels, the ground state |J), the 
resonance \R), and the final ionic state plus Auger electron \Fe) 



A{t) 



^ ai{t) ^ 



(33) 
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and the matrix 



H(t) 



DR{t)-'^W{t) ER-'^[VA + T,{t)]-u , (34) 

The final set of equations derived stepwise in the preceding sections now takes on the fol- 
lowing compact form 

iA{t) = Yl{t)7,{t). (35) 

The matrix H(i(:) can be viewed as the effective Hamiltonian governing the electron dynamics 
in the RA process of an atom in an intense laser pulse. 

Below we summarize the physical meaning of each term of the Hamiltonian matrix (13^ . 
The ground state energy Ej is augmented by the imaginary term —^Tph{t) describing the 
leakage from the ground state due to its direct photoionization (GS-leakage; see Eq. ( JT8l) ). 
Obviously, this term is time dependent because this leakage is only present during the pulse. 
Similarly, the total leakage from the resonance with dressed energy E^ — uj is provided by the 
imaginary part — |[r^ + r*(t)] on the respective diagonal element of H(t). The first term, 
— |r^, describes the usual time-independent leakage due to the Auger decay (the Auger 
decay rate or RA-leakage; see Eq. (jO])). The second term, — |r*(t), is the time-dependent 
leakage due to the direct photoionization of the resonance (RD- leakage; see Eq. (131]) ). 

Particularly interesting is the finding that the usual direct coupling D^it) between the 



ground state and the resonance through the laser field ([19|-|21|; see Eq. ( I5a|) ) is augmented 
by an additional a priori unexpected time-dependent term —^W{t). This term appears only 
if the photoionization of the ground state and the Auger decay are simultaneously treated 
(see Eq. fl2^ ). We have called this additional coupling Leakage-Induced Complex coupling 
(LIC-coupling). This coupling is an indirect coupling between the resonance and the ground 
state mediated by the combined action of the Auger decay and photoionization of the ground 
state and can be interpreted as follows: The photoelectron emitted by the ground state is 
recaptured by the residual ion to produce the resonance state (as described by the term 
— ^W{t) = —i7TD£,{t)V^ in Eq. fl24b|) ) and reversely the Auger electron can be captured 
by the residual ion which then becomes the neutral atom in its ground state (as described 
by the term -^W^{t) = -inVDl^it) in Eq. (Mah ). These processes mediate a coupling 
between the resonance and the ground state as long as the pulse is on. This coupling damps 
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the Rabi oscillations induced by the usual direct coupling DR{t) between the ground state 
and the resonance through the laser field. 

Finally, the last row of the Hamiltonian matrix (13^ describes how the studied specific 
final ionic state \Fe) is populated by the direct photoionization of the ground state (matrix 
element D£,.{t); see Eq. (11 II) ) and coherently by the Auger decay of the resonance (matrix 
element V^-; see Eq. (ISall ) all at a given kinetic energy ej of the emitted electron. The zeros 
in the last column of the Hamiltonian matrix make clear that there is no direct coupling 
between the ground state and resonance on the one hand side and the final ionic states on 
the other. One can first compute the dynamics of the coupled ground state and resonance 
and then subsequently determine the electron spectrum (see also Eq. fl24cp ). 



III. RESULTS AND DISCUSSION 



To exemplify the general theory of resonant Auger in intense X-ray fields we investigate 
the RA decay of the Ne Isi2s22p63p^ state (excitation energy is 867.12 ± 0.03 eV |28|). 
For this example all input data is available or could be computed, and the measurement 
is feasible. The Auger decay of this resonance in the XFEL-field has recently been studied 
theoretically in Ref. [l9|] where the isolated resonant channel populating the final ionic states 
has been considered. Here, we would like to study the impact of all possible mechanisms 
indicated in Fig. [T] We concentrate on two RA decay channels of different physical nature: 
the spectator RA decay into the ls^2s^2p'^(^D)3p and the participator decay into the 
ls^2s^2p^ final states. The corresponding features are clearly resolved in the RA spectra 
as their binding energies are 55.83 eV and 48.48 eV or, equivalently, the energies of the 



emitted electrons are around 811.29 eV and 818.64 eV, respectively 
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3Q|. 



The main difference between these two states consists in the very different ratios between 
the partial probabilities of the resonant and direct channels to populate them. This ratio 
is proportional to x -qty^- States populated by the spectator decay, briefly denoted 
as spectator states, are states with two missing electrons in the occupied and one excited 
electron in the unoccupied orbitals. They are typically seen as weak satellite lines in the 
direct photoionization spectrum (in XPS spectrum) because their production by a single 
photon is forbidden in the one-electron approximation and they appear due to electron 
correlation and relaxation effects. On the other hand, their partial Auger rate is usually 
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large since only occupied orbitals are involved in the Auger decay. Thus, for the spectator 
channel the ratio x is large, the Auger channel dominates over the respective direct ionization 
channel and the interference effects between these channels are expected to be small. In 
contrast, the participator states are main lines in the single photon XPS spectrum and, in 
addition, their partial Auger rates are small since it is the excited electron which is involved 
in the decay. As a result, the ratio x for a participator state is typically much smaller than 
that for the spectator state, and distinct interference patterns in the RA spectra can be 
expected. 

In order to illustrate the individual contributions of the different mechanisms to the RA 
spectra, the present calculations were performed not only employing the full formalism, but 
also employing systematically approximations which take into account only one or two mech- 
anisms. We shall address the results of the full calculations as 'exact' and the approximate 
ones as: 

• ^ Resonanf - only the resonant channel is accounted for as described in Sec. Ill A| 

• ^Direcf - only the direct ionization channel is taken into account as described in 
SeclllBl 

• ^Interference^ - both the direct and resonant channels and the interference between 
them are taken into account as described in Sec. Ill C\ 

• ^ExacV - all terms are taken into account including the direct photoionization of the 
resonance as described in Sec. IIIDI (see also Sec. Ill E|) . 

The quantities utilized in the present calculations are collected in Tab. [H The calculations 
are performed for the Gaussian-shaped pulse of a duration r centered at to 

(^(t) = e-(*-*o)'/^'. (36) 

In order to evaluate the LIC-coupling W{t), one has to include the manifold of all final 
ionic states (sum over index j in Eq. f l2^ ). As is clear from the above, the individual 
contributions of all states to W{t) are different in particular those of the spectator and 
participator states. Using the data in the literature, we were able to estimate only the 
modulus of the individual contributions of only one spectator and one participator state. 
In the absence of further data we assumed the same absolute values for all the important 
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TABLE I: Parameters for photoionization (PI) and resonant Auger (RA) decay of Ne utilized in 



the present calculations. The 


energy 


of the exciting 


; radiation is chosen to be = 867.12 eV. 


State/Quantity 


Value 


Source 


Method 


Ground state ls22s22pS 










Total direct PI cross section 


(7D = 


0.024 Mb 


T3„f [971 


Exerinient 


Resonance Isi2s22p63pi 










Excitation energy 


Er = 


867.12 eV 


Ref f28l 


T^i'y n PTi in pti t, 


Total RA decay rate 




270 meV 


Ref [281 


Experiment 


Oscillator strength 


Ir = 


0.0077 a.u. 


PypQpri-t- pcilr'nlo't'iri'nQ 


XXJ_ WiUii i t^iclAai biUii 


Total direct PI cross section 




0.025 Mb 




ij(7ii" n taI ti "vti "i" 1 on 


Participator state ls^2s^2p^ 










Energy of state 


Ef = 


48.48 eV 


Ref. [29] 


Experiment 


Branching ratio for RA decay 


X 


~ 2.1%") 


Ref. [29] 


Experiment 


Partial RA decay rate 


^ A 


— 5.7 meV 




Combination of hsted data 


Partial direct PI cross section 


^part 


= 0.015 Mb 


Ref. [27] 


Exeriment 



Spectator state ls22s22p'^(iD)3pi 

Energy of the state Ep — 55.83 eV 

Branching ratio for RA decay ^spct _ io.6%^) 
Partial RA decay rate 

Ratio of direct PI cross sections aspct/o-part='i% 

Partial direct PI cross section (T^pct = 0.0006 Mb 



pjct ^28.7 meV 



Refs. [29, 30] 
Refs. [29, 30] 

Ref. [311 



Experiment 
Experiment 

Combination of listed data 
Experiment (at 900 eV) 
Combination of listed data 



°) Present estimate from the Fig. 1 of Ref. 2^ . 



The intensities of the spectator transitions listed in Tab. 1 of Ref. 



29| are given relative to each other 



(i.e., they sum up to 100%) and not relative to the total Auger rate of 270 meV what is needed here. Since 



32% of the total intensity of the RA spectrum is associated with shakeup processes 

have been corrected accordingly. 



29( 1 ■ the values listed 



spectator and participator channels shown in Fig. 1 of Ref. [29| and have further neglected the 
possibility that different contributions could have different signs. The resulting value of the 
LIC-coupling can be considered to be an upper bound for the true value. Our computations 
shown below illustrate that for Ne the impact of the LIC-coupling is moderate for the field 
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intensities investigated. One can, however, anticipate that for other atoms and in particular 
for molecules the values and impact of the LIC-coupling can be substantial. 



A. Total electron yield 

The total electron yield as a function of the X-ray peak intensity calculated 'exactly' and 
compared with that obtained in the different approximations discussed above is depicted in 
Fig. [2] for the two pulse durations of r = 1 fs and 3 fs. The total electron yield was computed 
as in [l9| using 



lim / deAapAt^eAl'^ + lim / de^Aap, {t,e^^\'^ = 1 - lim |a/(t)p. (37) 

t-^-OO-^— ' J ■' ' ■' t^OO^—^ J ^ 3 ^ t^OO 



In the Resonant approximation (dash-dotted curves), only the decay of the resonance (RA- 
leakage) and the Rabi oscillations between the resonance and the ground state are competing 
with each other. This situation is discussed in [19]. If the duration of the pulse is comparable 
or shorter than the Auger decay lifetime of the resonance (r^ ~ 2.5 fs) and its intensity is 
low, the atom has a finite probability of staying neutral, i.e., of remaining in its ground 
state, after the laser pulse is over. As the intensity increases this probability decreases and 
the atom can be completely ionized with probability of 1 (see the first maximum in the 
total electron yield in Fig. [2]). At certain intensities the atom manages to complete several 
Rabi cycles during the pulse and the ionization probability drops again. These intensities 
correspond to the minima in the total yield 1^ , shown by the dash-dotted curves in Fig. [21 
and depend on the pulse duration (compare the upper and lower panels of this figure). 

The total electron yield computed in the Direct approximation (dotted curves) illustrates 
the individual contribution of the GS-leakage mechanism into the photoionization process 
of the atom. The X-ray pulse with energy below the Is-threshold may ionize the 2s and 
2p electrons of Ne by a single photon absorption. The total cross section for ionization of 
the valence shells {an = 0.024 Mb, Tab. [I]) is much smaller than the probability for the 
ionization via the resonant channel. Therefore, the total electron yield saturates to 1 due 
to the GS-leakage at much larger intensities than in the Resonant approximation. Because 
the ground state and the final ionic state posses different numbers of electrons, there are no 
oscillations in the total yield (note the absence of coupling between these states in Eqs. ( |T9l) ). 

In the Interference approximation (dashed curves in Fig. [21), three competitive mech- 
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FIG. 2: (Color online) Total electron yield after exposure of Ne to a coherent Gaussian-shaped 
pulse of duration 1 fs (upper panel) and 3 fs (lower panel). Shown are 'exact' results (solid 
line) and results of different approximations (broken lines) discussed in the text. These include the 
contribution of only the resonant channel {Resonant; dash-dotted line), of only the direct ionization 
channel {Direct; dotted line), and of both the direct and resonant channels and the interference 
between them {Interference; dashed line). The upward and downward arrows in the lower panel 
indicate the peak intensities chosen for the calculation of the Auger spectra. 
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anisms are present in the ionization of the atom. These are the decay of the resonance 
(RA-leakage) and the GS-leakage discussed above, and the interference between these two 
channels. As discussed in the theory, the interference is now modified from the weak field 
case by a complex term (LlC-coupling), as can be seen by inspecting the coupling in Sees. Ill CI 
and III El The whole coupling causes the Rabi oscillations of the population between the 
ground state and the resonance, and since it has now become complex, these oscillations 
can be modified as well and are usually suppressed. At those times when the population is 
in the resonant state (half-completed Rabi cycles) the RA-leakage mechanism is responsible 
for the ionization of the atom. On the other hand, if the population is in the ground states 
(completed Rabi cycles), the atom is ionized by the GS-leakage mechanism. The compe- 
tition between these two leakage mechanisms results in distinct modifications of the yield: 
the oscillations are now less pronounced and saturate along the trend imposed by the Direct 
approximation (cf the dashed-dotted and dashed curves in Fig. [2]). Consequently, the oscil- 
lations are suppressed from below. In addition, as can be seen in Fig. [21 the oscillations are 
also slightly suppressed from above and not all of them arrive at unity, and this weak but 
qualitatively interesting effect is a consequence of the complex coupling term (LIC-coupling) 
mentioned above. 

The RD-leakage mechanism enhances the ionization of the atom {Exact calculations, solid 
curves in Fig. |2]). Obviously, its impact is largest at the half-completed Rabi cycles, when 
the population is mainly in the resonance state, and is reduced by this leakage by ionization 
into other final ionic states thus suppressing somewhat the RA decay. The oscillations in 
the total electron yield now become even less pronounced and saturate much earlier as a 
function of the laser intensity (cf, dashed and solid curves). We note that depending on the 
pulse duration and the parameters of the system, there is a certain interval of intensities 
where both the GS-leakage and RD-leakage mechanisms are comparatively weak. In the 
present case of the longer pulse, these are the intensities up to the first minimum in the 
electron yield, before the atom manages to complete one Rabi cycle. 

B. Spectator Auger decay 

In order to illustrate the effect of the interference between the resonant and direct ioniza- 
tion channels we have computed the RA spectra for spectator and participator Auger decay 
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states. The calculations were performed for the full equations and also for all systematic 
approximations discussed above employing a pulse duration of r = 3 fs. The peak intensities 
chosen for the calculations correspond to the first three maxima and two minima in the total 
yield and also to the intensity related to the 50% of the yield as indicated by the upward 
and downwards arrows in the lower panel of Fig. [2J 

We first discuss the case of spectator RA decay. As was mentioned above, the ratio of the 
partial probabilities for the resonant and direct channels, x d^Va ' ^^^S^ for spectator 
RA decay. Therefore, one intuitively expects here only small or moderate interference effects 
between the resonant and direct photoionization channels. The RA spectra computed for 
the spectator channel are depicted in Figs. [3] and ID As was shown in Ref. 19|, the Rabi 
flopping induced by a coherent X-ray pulse leads to a strong modification of the purely 
resonant Auger electron line profile (dash-dotted curves in the figures). With the increase of 
the peak intensity, the single RA electron line found at rather weak fields bifurcates resulting 
in local minima in the line shape. This effect was also discussed in a different context 25 
years ago [32 1. 

The individual contributions of the partial direct photoionization of the ground state to 
the spectral line discussed are shown by the dotted curves in the figures (note that these 
curves are shown on an enhanced scale compared to the other curves indicated by a factor 
xK). Since there is no coupling between the ground and the final states (see Eqs. f|T9|) ). 
the line profile computed in the Direct approximation remains unchanged with increasing 
laser intensity and hence displays no bifurcations. The intensity of the line in the spectrum 
increases rapidly as the field intensifies. Surprisingly, at peak intensities around the second 
minimum and the third maximum in the total electron yield (see Fig. |2]), the individual 
contribution of the direct channel becomes about the same order of magnitude as that 
of the resonant channel as seen by comparing the dash-dotted and dotted curves in the 
uppermost spectra of Figs. |3] and HI This finding can be rationalized. As seen from Fig. |2l 
the effect of the resonant channel rises fast with the intensity of the field and except of the 
impact of the Rabi oscillations saturates already at the first maximum of the total electron 
yield. The effect of the direct channel, on the other hand, becomes relevant only at higher 
field intensities and saturates much later. As soon as the latter is also saturated, the effect 
of the direct channel on the spectrum becomes comparable to that of the purely resonant 
channel. 
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FIG. 3: (Color online) The spectator Auger electron spectra computed for a pulse duration of 3 fs 
employing the full equations ('exact') and the systematic approximations discussed in the text. 
The peak intensities chosen for the calculations are indicated by the downward arrows in the lower 
panel of Fig. [2j Note that the dotted curves are shown on an enhanced scale compared to the other 
curves indicated by the factor xif at the right hand side of the curves. 

The interference between the resonant and direct channels (Interference approximation; 
dashed curves) results in a slight asymmetry of the RA line profiles for the spectator state. 
With the present parameters for the specific transition studied here, the interference is de- 
structive on the low electron energy side, and constructive on the high energy side. Whether 
the interference is destructive or constructive depends, of course, on the relative signs of the 
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FIG. 4: (Color online) The spectator Auger electron spectra computed for a pulse duration of 3 fs 
employing the full equations ('exact') and the systematic approximations discussed in the text. 
The peak intensities chosen for the calculations are indicated by the upward arrows in the lower 
panel of Fig. [2j Note that the dotted curves are shown on an enhanced scale compared to the other 
curves indicated by the factor xif at the right hand side of the curves. 

corresponding amplitudes. In the present model we assumed that the resonant and the di- 
rect amplitudes have the same signs. If these signs were different, the computed spectrum 
would be just the mirror reflection of the line profile at the center. 

Obviously, the role of the direct photoionization, and, hence, of the interference effects 
increases with the field strength. However, as the total electron yield saturates as a conse- 
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quence of the GS-leakage mechanism, also the impact of the interference effects saturates. 
As the field increases further, the spectra become overall suppressed due to the depopula- 
tion of the resonance into other usually high lying final ionic states (leakage of the resonance 
state). The impact of this RD-leakage is visible in the resonant Auger spectra as can be 
seen in the figures by comparing the results of the full calculations (solid lines) with those 
of the interference approximation (dashed lines) not including this leakage mechanism. 



C. Participator Auger decay 

The 'exact' and approximate theoretical RA spectra of the participator channel are de- 
picted in Figs. |5] and |6l The spectra computed in the Resonant approximation show the 
same trends as for the spectator channel (dash-dotted curves). Because of the smaller par- 
ticipator Auger rate (5.7 meV compared to 28.7 meV for the spectator channel; see Tab. [11), 
the participator RA electron peaks are by about 5 times weaker than those found in the 
spectator spectra at the same peak field intensities. 

At the same time the partial direct photoionization cross section of the participator 



channel is 25 times larger than that of the spectator channel [Sjj {(Jpart = 0.015 Mb and 
o'spct = 0.0006 Mb; see Tab. [T]). Consequently, the individual contribution of the direct 
channel to the RA peaks is now seen already at relatively low field intensities {Direct ap- 
proximation; dotted curves in Figs. O and E]). The dotted curves in these figures are decreased 
by a factor 1/K compared to the other curves in the figures. The electron spectra computed 
in the Resonant and Direct approximations are already comparable in magnitude at the 
intensity at which 50% of the total ion yield is achieved for the first time in Fig. [2l This is 
clearly seen when comparing the dash-dotted and dotted curves in the lowermost spectrum 
in Fig. O Obviously, the role of the direct channel increases greatly with the field strength as 
was also the case in the spectator channel. If the direct photoionization and resonant spectra 
were independent of each other, the direct ionization would be the dominant mechanism for 
the population of the participator final ionic state at the lager field intensities shown in the 
figures (note the factor l/K). 

However, interference effects are found to be substantial and the respective spectra ex- 
hibit much interesting and distinct structure. The RA spectra computed in the Interference 
approximation are depicted in Figs. and E] where they are shown as by dashed curves. By 
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FIG. 5: (Color online) The participator Auger electron spectra computed for a pulse duration of 
3 fs employing the full equations ('exact') and the systematic approximations discussed in the text. 
The peak intensities chosen for the calculations are indicated by the downward arrows in the lower 
panel of Fig. [2l Note that the dotted curves are shown on a suppressed scale compared to the other 
curves indicated by the factor 1/K at the right hand side of the curves. 

comparing these curves with the purely resonant spectra (dash-dotted curves), one imme- 
diately concludes that the interference between the resonant and direct pathways for the 
population of the final participator state is very important at all the studied intensities of 
the field. The interference patterns associated with the asymmetry of the RA spectrum are 
strongly pronounced. As the consequence of the present parametrization of the amplitude's 
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FIG. 6: (Color online) The participator Auger electron spectra computed for a pulse duration of 
3 fs employing the full equations ('exact') and the systematic approximations discussed in the text. 
The peak intensities chosen for the calculations are indicated by the upward arrows in the lower 
panel of Fig. [2l Note that the dotted curves are shown on a suppressed scale compared to the other 
curves indicated by the factor 1/K at the right hand side of the curves. 

signs, one observes destructive and constructive interferences on the low and high electron 
energy sides of the spectra, respectively (see also discussion in the preceding subsection). 
We also notice here that the fingerprints of the interference in the ls^2s^2p® participator 
channel might be visible even in the weak field RA spectrum measured with conventional 
synchrotron radiation: the slight asymmetry of the peak at binding energy of 48.48 eV in 
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Fig. 1 of Ref. 291] is very similar to that in the lowermost spectrum in Fig. 5. 

Moreover, the interference determines not only the shape of the computed RA spectra, 
but their integral intensities as well. Interestingly, the integral line intensity computed in 
the Interference approximation at the larger intensities of the field studied (middle and 
uppermost spectra in Figs. \5\ and E]) is larger than that in the Resonant approximation, but 
much smaller than that in the Direct approximation (note factors | and ^) . This is due 
to the competition between the Auger decay and GS-leakage mechanisms and persists also 
in the 'exact' calculations where the leakage from the resonance is taken into account (solid 
curves). Finally, we note that the impact of the latter losses to higher ionic states (RD- 
leakage mechanism) on the decay spectra is to suppress their intensity without substantially 
changing their line shape. 



IV. CONCLUSIONS 



The theory of resonant Auger decay of an atom exposed to a coherent monochromatic 
X-ray pulse of strong intensity is presented. The theory incorporates the resonant and direct 
photoionization pathways from the ground electronic state to the final ionic states, as well 
as the ionization losses of the resonance state. Interesting interference effects appear in the 
electron spectra due to the presence of the direct ionization pathway and the excitation and 
decay of the resonance which populate the final ionic state coherently. Due to the leakage 
of the ground state, the coupling between the ground electronic state and the 'dressed' reso- 
nance state is modified giving rise to an additional complex coupling term which we call the 
Leakage-Induced Complex coupling (LIC-coupling). All together, the effective Hamiltonian 
governing the production of ions and electrons is non-hermitian where all its terms, diagonal 
and off-diagonal, are complex. 

The theory predicts a strong competition of several mechanisms of ionization of the atom. 
While the pulse proceeds, the non-hermitian coupling induces (damped) Rabi oscillations of 
the population between the ground state and the resonance. While the population is in the 
ground state, in particular when a Rabi cycle is completed, it is subject to leakage into ionic 
states via direct photoionization (GS-leakage). In the time intervals when the population is 
in the resonant state, in particular at half-completed Rabi cycles, it can decay by resonant 
Auger decay (RA- leakage), but it can also be transferred into highly-excited ionic states via 
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direct photoionization of the resonance itself (RD-leakage). 

The developed formalism is applied to study the excitation and Auger decay of the 
ls^2s^2p^3p^ resonance of a Ne atom in the intense field of an X-ray pulse with a central 
frequency of 867.12 eV. It is shown that at relatively weak intensities of the field, the resonant 
Auger decay is the dominant mechanism of ionization of Ne, but as soon as this channel of 
ionization saturates with the increase of the field intensity, the GS-leakage as well as the 
RD-leakage mechanisms start to be competitive. At larger intensities, when the impact of 
the GS-leakage and RD-leakage saturates as well, all these three different mechanisms of 
production of ionized Ne become of comparable magnitude. 

The interference effects are of different importance to spectator and participator final 
states of the Auger decay. To demonstrate the interference effects between the resonant and 
direct photoionization pathways, we have studied the electron spectra for the two specific 
cases of the spectator ls^2s^2p^(^D)3p and the participator ls^2s^2p^ final Auger 
states. The interference effects increase with the field strength until they saturate when 
the total electron yield saturates due to the GS-leakage mechanism. For the final state 
produced by the spectator decay, the decay by resonant Auger is the dominant ionization 
mechanism, and the interference results in a slight asymmetry in the line profile of the 
spectrum. In contrast, the direct photoionization is the dominant mechanism at large field 
intensities in the case of participator final states. Nevertheless, a distinct interference pattern 
appears in the electron spectrum which remains highly structured in spite of the strong direct 
photoionization at high field intensities. This finding makes the study of resonant Auger 
spectra in strong fields worth while. 

We would like to conclude with the following remark. In molecules, the discussed com- 
petition of different leakage mechanisms and the interference effects will be much more 
intricate, due to the presence of the nuclear dynamics. In particular, we expect in molecules 
also the appearance of strong non-adiabatic effects induced by the field [33|, l3J]- Here, we 
suggest to follow the emerging competition and interference effects in the time domain which 
will hopefully be possible with the advances in the development of pump-probe techniques 
utilizing femtosecond and sub-femtosecond pulses. 
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Appendix A: Local approximation 

In this appendix we discuss the local approximation which leads to Tph{t) in Eq. f[T5]) as 
an example. The starting point is the set of equations f|T3|) . The straightforward integration 
of this set is impossible owing to the continuous set of Eqs. ( 113b|) describing the final ionic 



state together with the photoelectron. In order to evaluate the integral in Eq. fll3ap . one 
would have to solve the whole set of Eqs. (]13b|) simultaneously for all electron energies which 
is a formidable task. The problem becomes tractable by implying the local approximation 
, [25 1 to integrate Eqs. (11 5p and f ll6p . Due to the non-locality of the kernel fll6p . the 
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amplitude ai{t) is needed at all times smaller than t, i.e., the kernel is a so-called 'memory 



kernel' 



25j. The local approximation consists in eliminating this 'memory' by loca. 



ization 



of the kernel (fT6|) . In order to proceed, we make the same assumptions as in Ref. 26| where 
the local approximation is derived in detail. 

At first, we assume that the integral function D\j{t)Dr)it') in Eq. (|T6l) varies negligibly as 
a function of energy in the symmetric energy interval of (^o ~ ^0 + e) around the position 
of the center of the spectral line of the emitted electron sq = Ej + u — Ep, and vanishes 
otherwise (the choice of the interval length, parameter e, will be justified below). The 
approximation is good as long as Eq is large, i.e., the ejected electron has a large energy and 
the ionization process is far from the ionization threshold. In this approximation we obtain 



K{t -t')= I D^^it) DB(t')e-'(^^+"°-")(*-*')e-'(^-^°)(*-*')rf£ = 
D^j,{t) Doit') e-^(^^+-o-c.)(t-t') r ^~Mt-t')^ 



2sin|e(t^^ tO] ^t^(^)^^^^.)^-.(E,+.o-^)(t-t')_ (Al) 

Next, we assume that the time-dependence of the amplitude ai(t) can be factorized in two 
parts of different physical origins 

ai{t) =e-'^''af\t). (A2) 
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The phase factor e~*^^* describes the time- dependence of the amphtude owing to the energy 
of the ground state. The rest of the time-dependence of the amphtude a/(t) is associated 
with the impact of the laser pulse, and is deposited in the function af\t). 
The time integral in Eq. f[T5]) can now be simplified 

/t /»oo 
K{t-t')ai{t')dt' = -i / K{T)ai{t-T)dT = 
-oo Jo 

Jo ^ 

-zDUt)e-''''' r ^^^Doit - r)af\t - T)dT, (A3) 
Jo ^ 

where in the last step we used the definition of the center of the electron spectral line 
Ef+Eq — u — Ej = 0. The main difference in further integrating Eq. ( 1A3I) from the procedure 



described in Ref. 



26| consists in the time dependence of the DD(^)-function. To proceed 



we note that the timescales at which the functions in the integrand undergo substantial 
variations depend on the pulse duration and on the yet to be determined parameter e. The 
function D£){t — T)af\t — T) in Eq. f lASp varies on the timescale of the pulse duration (usually 
a few femtoseconds). The relevant timescale for substantial changes of the other appearing 
function, i^l^tll^ is determined by the parameter e and is of the order of In order to be 
able to perform the integration of Eq. ( IA3I1 over time analytically, one has to make the latter 
timescale much shorter than the duration of the pulse. For this purpose, the parameter e 
has to be taken as large as possible. 

The upper limit for the parameter e is naturally determined by the ionization threshold 
(e ^ ^o)- III the case of core-excitation, the electron energy Eq is large (in the present case 
more than 800 eV). By taking the parameter e = Eq (as in Ref. 26|), the timescale for 
substantial variations of the function ^'^'"^^^^ can be made much shorter than the duration 

r 

of the pulse. For the present case, this timescale is in the sub- femtoseconds domain (|^ ^ 
0.005 fs). Finally, we note that the function ^'^"^^'^"^^ vanishes rapidly as r increases. As a 



result, the main contribution in the integral flASp comes from the times around t ^ 0. The 
analytical integration of Eq. (IA3I1 is now straightforward and we obtain 



- iDUt)DD{t)e-'^''af\t) [ ^^^^^dr = -Z7r|Z}z)(t)| V(t), (A4) 

Jo ^ 
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and hence Eq. ( !T7|) as the final result. 
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